Detecting Features in the Dark Energy Equation of State: A Wavelet Approach 
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We study the utility of wavelets for detecting the redshift evolution of the dark energy equation 
of state w{z) from the combination of supernovae (SNe), CMB and BAO data. We show that local 
features in w, such as bumps, can be detected efficiently using wavelets. To demonstrate, we first 
generate a mock supernovae data sample for a SNAP-like survey with a bump feature in w{z) hidden 
in, then successfully discover it by performing a blind wavelet analysis. We also apply our method 
to analyze the recently released "Constitution" SNe data, combined with WMAP and BAO from 
SDSS, and find weak hints of dark energy dynamics. Namely, we find that models with w{z) < — 1 
for 0.2 < z < 0.5, and w{z) > —1 for 0.5 < z < 1, are mildly favored at 95% confidence level. This 
is in good agreement with several recent studies using other methods, such as redshift binning with 
principal component analysis (PCA) (e.g. Zhao and Zhang, arXiv: 0908.1568). 



I. INTRODUCTION 

The observed acceleration of cosmic expansion [1, 2] 
means that either our understanding of gravity on cos- 
mological scales is incomplete, or that there exists a yet 
unknown dark energy (DE) component, which is gravita- 
tionally repulsive, contributing around 70% of the total 
energy budget of the universe. The redshift evolution of 
the DE fraction is determined by its equation of state 
w{z) via 

nj,Eiz) = nuECxp{-3 J dlnz'[l + w{z')]} , (1) t 

where f^DE is the dark energy fraction at z = 0. The 
simplest DE candidate is the cosmological constant A 
with w = —1. However, in the absence of a convincing 
theory of DE, one should not disregard the possibility of 
a more complicated evolution history of w{z). 

Much effort has been put in extracting w{z) from data. 
Even though A provides a good fit [3] , current experimen- 
tal constraints on w are not accurate enough to rule out 
the possibility of a time-varying w. It is of fundamental 
importance to know if DE is a constant or dynamical, 
and future observations will be in a much better posi- 
tion to address this important question [4]. Approaches 
used so far to estimate w{z) from data can be classi- 
fied into three types: fitting a model [5-8], binning and 
principal component analysis (PCA) [9-11], and a direct 
differentiation of distance data [12]. In the first method, 
one assumes a functional form for w(z) or the Hubble 
parameter H{z) and fits the form to data. Many differ- 
ent functions have been considered and the results are 
generally strongly dependent on the assumed model. In 
the PCA method, recently adopted by the JDEM Figure 
of Merit Science Working Group [13], wlz) is binned in 
redshift or scale-factor, and the value in each bin is con- 
sidered as a parameter of the model. The best measured 
principal components can be used to reconstruct w{z). In 
the third method, the Type la SNe luminosity distance 
data is smoothed, and w{z) is calculated from the second 
derivative of this smoothed function. 



All of the three techniques above have their merits, and 
can be used to address appropriately formulated ques- 
tions. In this paper we are examining the possibility of 
using wavelets as a tool to detect possible features in 
w{z) from a combination of supernovae (SNe), cosmic 
microwave background (CMB) and baryon acoustic oscil- 
lations (BAO) data, in a model-independent way. Local 
featues in w{z) would exist if, for example, dark energy 
was a quintessence scalar field with a potential that had 
a a step-like bump. Sharp features in w{z) are also pre- 
dicted in certain types of phantom models of dark energy 
14]. 

An advantage of wavelets is that one does not need 
to assume a prior knowledge of the position and size of 
features in w(z). Wavelets search for large and small scale 
features at the same time and can reconstruct them with 
only a few parameters. 

This paper is organized as follows. In Sec. II we de- 
scribe the wavelet method. In Sec. Ill we demonstrate the 
method by hiding a feature in a simulated data and then 
performing a blind analysis to reconstruct it. In Sec. IV 
we apply the same method to search for features in the 
current data and conclude with a summary in Sec. V. 



II. METHODOLOGY 
A. The wavelet parametrization 

Wavelets, known as "mathematical microscopes" , have 
the ability to "see" signals on multiple scales, making 
them a powerful tool for detecting features in noisy data. 
Examples of their use in cosmology include reconstruct- 
ing the primordial power spectrum [15, 16], probing the 
non-Gaussianity of the primordial fluctuations [17], and 
searching for signals of parity breaking [18]. Wavelets are 
specially designed functions that have the property of be- 
ing localized in both frequency and configuration space. 
Individual wavelets are generated from a "mother" func- 
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tion as 



(2) 



where is the mother function, L is the range on 

which the function is defined, < m < 2" — 1 and 
< n < — 1. Each wavelet function has a locahzed 
shape. The index n determines the size of the feature, 
while m determines its location on x. Functions ipnm arc 
orthogonal by design, hence we can expand any function 
onto the wavelet basis in the same way as we do for any 
other orthogonal basis: 



2"-l 
n— m— 



(3) 



By this definition, each wavelet acts as a band-pass filter 
and an infinite number of wavelets are needed to cover the 
whole spectrum. This problem is solved by introducing a 
scaling function. Wavelets coefficients together with the 
coefficient of the scaling function, 2^ in total, can be 
used to describe f{x) at 2^ points. 




FIG. 1: The 5th and 15th D4 (upper) and Haar (lower) 
wavelets. The resolution (16 points) in this illustration is the 
same as the one used in the numerical calculations in this 
paper. 

In this work, we consider w{z) to be an unknown func- 
tion which can, in principle, have local features in z. Like 
any other function, we can expand it into the wavelet 
basis as 



l+w{Zj) = ^P,'lp,{Zj) 



(4) 



where Zj are the redshift points at which we calculate 
If, Pi are the wavelet expansion coefficients, and ipi are 
the wavelet functions. We chose to simplify our indices 
so that Pi is the scaling function coefficient, P2 = oooj 
Pa = fliOj etc. We then treat the wavelet expansion co- 
efficients as the parameters of the model to be measured 
from data. Note that we expand 1 -I- u'(z), as opposed to 
w{z), so that all coefficients are zero if w = —1, as in 
the currently favored ACDM model. We use Daubechics 



4 (D4) [19] and Haar [20] wavelets (prototypes shown in 
Fig. 1) in this paper. The D4 wavelets are highly localized 
and suitable for reconstructing sharp features. The Haar 
wavelets are made of step functions. We should keep in 
mind that our results using different types of wavelets 
should converge and be consistent with each other in the 
limit of large number of data points. 

Based on the expected redshift range of a future Type 
la SNe survey like SNAP, we take the redshift range to be 
< 2 < 1.7. We choose to work with up to 16 wavelet pa- 
rameters, which allows us to calculate w{z) at 16 equally 
spaced points in z and probe features of up to 1.7/8 ~ 0.2 
size in redshift. The values of the wavelet coefficients 
are determined by fitting (4) to data using Markov-chain 
Monte Carlo (MCMC) as explained in Subsection II B. 



B. How many wavelet coefficients? 

The main advantage of wavelets is that we do not 
need to a prior knowledge of the position, the scale and 
the shape of the possible features in the function w{z). 
Wavelets search for large and small scale features at the 
same time and the fact that they are localized in con- 
figuration space makes it possible to reconstruct features 
with only a few nonzero coefficients. In practice, we start 
with performing an initial MCMC run with all 16 wavelet 
coefficients as free parameters. We then keep those pa- 
rameters whose signal-to-noise ratio (SNR) is above a 
certain threshold, and set parameters with low SNR to 
zero, since they are not sensitive to features in w{z). From 
the surviving set of wavelet coefficients, we find the com- 
bination that has the lowest reduced x^. This gives us 
the optimal number of wavelets we need to detect a fea- 
ture in w{z). We then use three model selection criteria, 
described below, to determine if this wavelet fit is in fact 
preferred over a constant w. 

To judge whether a model is preferred by data, we use 
three selection criteria: the improvement in per extra 
degree of freedom Ax^/A(DoF), the Akaike Information 
Criterion (AIC) [21], and the Bayesian Information Cri- 
terion (BIC) [22]. The AIC and BIC are defined as 



AIC 
BIC 



-2\ogC + 2N„ 



2Np{Np + 1) 
Nd-Np-1 



-2 log/: 



NpNdlog{Nd) 
Nd-Np-1 



(5) 
(6) 



where Np is the number of parameters and iV^ is the num- 
ber of data points. The second terms on the right hand 
sides of these relations show the penalty for adding new 
parameters, which is stronger for BIC. Both expressions 
include second-order corrections, which is important for 
the case with a small data sample size, i.e. small Nd- The 
AIC and BIC are widely used in the model selection - the 
lower the value of AIC(BIC), the more the model is pre- 
ferred. A model with ABIC < —10 is decisively preferred, 
while models with -10 < ABIC < -6, -6 < ABIC < 
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z 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 


JDEM 


Niz) 

am{z) (xlO-^) 


300 35 64 95 124 150 171 183 179 170 155 142 130 119 107 94 80 
9 25 19 16 14 14 14 14 15 16 17 18 20 21 22 23 26 



TABLE I: The redshift distribution of type la supernovae N{z) for a SNAP-like JDEM survey [23], together with 300 SNe from 
the NSNF [24]. The redshifts given are the upper hmits of each bin. Magnitude errors am{z) are evaluated at bin midpoints. 



—2 and ABIC > 2 have very strong, strong, positive and 
weak evidence of model preference, respectively. In the 
case of the AIC, models with AAIC < 2 have substantial 
support, models with 4 < AAIC < 7 have considerably 
less support and those with AAIC > 10 have essentially 
no support. 

C. Current and future data sets 



Sm being the expected uncertainty and Zmax the maxi- 
mum redshift. Table I shows the redshift distribution of 
type la supernovae expected to be detected by JDEM 
together with estimated error in magnitude for different 
redshift bins. In addition, we assume that 300 supernovae 
at z < 0.1 will be measured by the Nearby Supernovae 
Factory (NSNF) [24]. For current SNe data, we use the 
recently released "Constitution" sample with 397 data 
points [25]. 



For our tests we will work with a combination of SNe, 
CMB and BAO data, considering both the current avail- 
able data and the simulated future data. The combined 
likelihood of the model is calculated as 



- 2 log £ = xtotal = XsNc + XCMB + Xbao (7) 



with 



X 



(x-d)^C-i(x-d) 



(8) 



where x is theoretical value for an observable, d is data 
vector and C is the data covariance matrix. 



1. SNe data 

The observable quantity in SNe data is the rcdshift- 
dependent magnitude: 



m = M + 5\ogdL + 25 



(9) 



where M is the intrinsic supernova magnitude and (11 the 
luminosity distance. The luminosity distance is defined as 



dL{z) = {l + z) 



^ dz' 



' 



(10) 



where H(z) is the Hubble parameter with a current value 
of Hq. We use the characteristics of a SNAP-like JDEM 
survey [23] to simulate futuristic SNe data with errors 
given by 



0-m{z) 



dm^ 



(11) 



where an 



0.15. The systematic error, dm, is assumed 



to increase linearly with redshift: 
dm — dm 



Zn 



(12) 



2. CMB and BAO data 

Rather than working with the full CMB data, we use 
the procedure suggested in [3] and work with two distance 
ratios: the shift parameter R, and the acoustic scale £aj 
defined as 



i?(z,) = {1 + z,)Da{z.) (13) 



(l + z*) 



TsiZf) 



(14) 



where z^, is the redshift of decoupling, Da{z) is the proper 
angular diameter distance to redshift z, and rs{z) is the 
comoving sound horizon at z. It was shown that these two 
quantities contain almost all of the information required 
to constrain dark energy with CMB [26-28]. For future 
CMB data, we used the Planck characteristics for param- 
eters and their errors [29]. For today's data we used the 
WMAP-5 values and errors for the parameters [3]. 

We use the BAO data point at z = 0.35 from [30] 
corresponding to r^/ Dy(z), where Dy{z) is given by 



Dv = 



{l + zfD\{z) 



H{z) 



1/3 



(15) 



where c is the speed of light. 



D. PCA analysis 

A significant departure from zero of any of the wavelet 
coefficients would indicate a departure from the ACDM 
model. Although the wavelet functions are orthogonal 
by construction, and hence their expansion coefficients 
should be uncorrelated, the error bars on the fitted 
wavelet coefficients are actually correlated. In the pres- 
ence of correlation among parameters, one can not deter- 
mine conclusively at what significance level a deviation 
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from zero has occured for any one of the coefheients. In 
order to quantify the detection of a departure from zero 
we decorrelate the wavelet coefficients. The PCA analysis 
[31] is designed to do precisely that by diagonalizing the 
covariance matrix of the original parameters to find their 
uncorrelated linear combinations. The original data vec- 
tor Pj is transformed into a new data vector qi using the 
same rotation matrix coij that diagonalizes the covariance 
matrix. Namely, we have 

q, = ^uj,jPj (16) 

3 

SO that the new covariance matrix, Cq = uj Cp uj^^ , is 
diagonal, and qiS are the uncorrelated principal compo- 
nents of P. Like the wavelet coefficients, principal com- 
ponents should also be zero for the ACDM model. We 
will use the PCA to quantify the significance with which 
the current data sees departures from ACDM. 

E. Features in w{z) 

Broad features in w{z), such as a single transition from 
one constant value to another, can be efficiently detected 
by other methods such as binning. Our aim is to de- 
velop a method to detect, if any, local features in w(z) 
in a model-independent way. In this sense, we found that 
wavelets can do a good job phenomenologically. Wavelet 
method is more effective for searching for localized fea- 
tures in w{z) that would likely require a larger number 
of parameters to describe with other methods. From the 
theoretical point of view, a local feature in w{z) is not 
particularly motivated, but can be produced in scalar 
field models of dark energy with a step-like bump in the 
potential, or in certain types of phantom models [14] . 

In the next section we will try to recover a feature in 
w{z), a bump, that we "hide" in a simulated future data. 
Such a bump is not ruled out by current data [32] , since 
the averaged 'w{z) - the quantity that the current data es- 
sentially constrains - does not significantly deviate from 
— 1 for this feature. The main rationale for considering 
this particular choice is that such a bump feature cannot 
be easily discovered by binning if one does not know the 
size, the position and the scale of the bump in advance. 

III. RECONSTRUCTING FEATURES WITH 
FUTURE DATA 

In this section, we use the wavelet method to recon- 
struct a feature hidden in the simulated mock future 
data. Our procedure can be summarized as follows: 

1. Generate a mock dataset with a feature of 'w{z) 
hidden in; 

2. Perform the first MCMC fit for 16 wavelet coeffi- 
cients; 
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FIG. 2: Comparison of the residuals of the distance modu- 
lus for ACDM, WCDM and bump models with characteristic 
error-bars for SNAP. 

3. Sort the wavelet coefficients according to their 
signal-to-noisc ratio (SNR) and set all except the 
few with the highest SNR to zero; 

4. Perform the second MCMC with all possible com- 
binations of "surviving" wavelet coefficients to find 
the model with the least per degree of freedom 
(DoF) 

5. Quantify the detection using Ax^/A(DoF), AAIC 
and ABIC. 

To start, we simulate future SN data for a dark energy 
model with a feature in u;(z) given by 

= -l+^exp[-(z-zo)'M , (17) 

where we choose A = —0.9, zq = 0.6, a = 0.1. This w{z) 
is shown with a black solid line in Fig. 5. We assume a fiat 
universe, and choose the other cosmological parameters 
as rjfe = 0.046, U.m = 0.279 and Hq = 70 km/s/Mpc. We 
also apply the distance priors derived from the CMB and 
BAO data as described in Sec. II C. We assume a SNAP- 
like JDEM survey with the luminosity data binned into 
17 redshift bins in the range of 2 G [0, 1.7]. Fig. 2 shows 
the comparison of the residuals (defined with respect to 
the ACDM model) of the distance modulus for the model 
with a bump-like feature in w{z) and a constant-w model 
which has the same average value w (i.e. one that gives 
the same distance to the CMB last scattering surface). 
Error-bars are plotted according to the characteristics of 
SNAP. The plot shows that a SNAP-like experiment is 
capable of differentiating between these models around 
z = 0.5. 

Having produced our mock datasets, we perform a 
MCMC fit for all the 16 wavelet parameters plus flbh^, 
flrnh'^ and h, assuming a fiat universe and flat priors for 
all the wavelet parameters. Fig. 3 shows the marginal- 
ized 1-D likelihood distributions for the 16 wavelet pa- 
rameters, with Table II listing the corresponding mean 
values and the 68% CL error bars. As one can see, for 
many wavelet parameters the standard deviation is much 
larger than the mean value, indicating that they are not 
sensitive to the presence of a feature in w. Such poorly 
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FIG. 3: The 1-D posterior distributions of all the 16 wavelet coefRcients obtained from our first MCMC run. 



constrained parameters are thus redundant, so we keep 
only the ones with the highest SNR. In order for the D4 
wavelets to be able to describe a constant w, one should 
keep Pi. In addition, we keep the 2-3 coefficients with 
the highest SNR. In our example, we keep parameters 
PI, P5 and P12. 

Next, we redo the MCMC using every combination of 
the selected parameters (with PI always included) and 
find the one that gives the minimum value of the reduced 
X^. In this case, the optimal combination is that of PI 
and P5, while setting P12 = 0. Marginalized likelihood 
distributions for this model are shown in Fig. 4 and the 
plot of w{z) reconstructed using this model is shown in 
the left panel of Fig. 5. 



Para. 


Constraints 


SNR 


Para. 


Constraints 


SNR 


PI 


-0.1 ± 1.2 


0.1 


P9 


-0.5 ±3.0 


0.2 


P2 


1.1 ±4.4 


0.3 


PIO 


-0.7 ±4.3 


0.2 


P3 


-1.3 ±4.3 


0.3 


Pll 


1.5 ±4.6 


0.3 


P4 


-0.2 ± 1.4 


0.1 


P12 


7.4 ±9.2 


as 


P5 


1.7 ±2.7 


06 


P13 


-0.1 ±5.0 


0.0 


P6 


-0.3 ±4.7 


0.1 


Pll 


0.1 ± 5.2 


0.0 


P7 


0.2 ±5.2 


0.0 


P15 


1.6 ±4.8 


0.3 


P8 


-0.3 ± 1.7 


0.2 


P16 


-1.6 ±4.1 


0.4 



TABLE II: Marginalized mean values with 68% CL. errors 
and SNR for all the 16 wavelet parameters. The high SNR 
parameters are underlined. 




0.0224 0,0228 0.0232 0.69 0.70 0.71 0.72 0.71 0.72 0.73 0,74 

n h £2. 



FIG. 4: 1-D posterior distributions of the 1st and the 5th 
wavelet coefficients, and the other four cosmological parame- 
ters. 



Model 



Ax /A(DoF) AAIC ABIC 



WCDM(w = -0.9) 
CDM-D4-P1,5 
CDM-Haar-P2,4 









-8.8 


-6.3 -3.7 


-8.3 


-6.8 -3.2 



TABLE III: The differences in the model selection criteria for 
the simulated data with a bump feature in w(z). 



In the above example two wavelet parameters were 
measured to be non-zero. Their uncorrclated linear com- 
binations, or principal components, have essentially the 
same values: 

PCI = 0.9 ±0.4 , PC2 = 0.43 ± 0.08 , (18) 



indicating evidence for ~6ct and ^2a departures from 
zero respectively. The correlation between the two 
wavelet parameters is negligible, which need not happen 
in general, but is somewhat expected. Namely, this com- 
bination is the one most favored by the reduced x^, which 
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rcdshift z 

FIG. 5: The hidden feature (black solid), and the reconstruc- 
tion using D4 (red dashed, left panel) and Haar (red dashed, 
right panel) wavelets. The horizontal dashed line and the 
shaded band illustrate the fit with the constant w model and 
the 68%CL errors. 



means that both parameters are likely needed to describe 
the data - one of them cannot compensate for the other. 
This implies a small correlation between these two pa- 
rameters. 

Detecting a departure from w = —\, which is all we 
have discussed so far, does not in itself amount to de- 
tecting the bump feature. The value of w{z) in Fig. 5 
averaged over redshift is less than ~1, so we would de- 
tect a departure from —1 by simply fitting a constant 
w. Hence, the more relevant question is not to what ex- 
tent wavelets are able to detect a departure from —1, but 
to what extent the wavelet fit is preferred to a fit by a 
simpler model with one constant w parameter. 

To answer this question, we fit the model with a con- 
stant w (WCDM) in a fiat universe and apply the model 
selection criteria to determine the degree to which the 
wavelet model is preferred in comparison. In particular, 
we calculate the change in per additional degree of 
freedom, the difference in the BIC and AIC. As one sec 
from the numbers in Table III, the wavelet fit is strongly 
preferred according to all three criteria. 

In addition to using D4 wavelets, we have followed the 
same procedure and tried to recover the bump feature 
using the Haar wavelets. We find again that a model with 
two wavelet parameters is optimal and is preferred to 
WCDM at a similar significance level as in the D4 case. 
The plot of 'w{z) reconstructed using two Haar wavelets 
is plotted in the right panel of Fig. 5. 

As a test, we have repeated the same procedure using 
data simulated according to parameters of existing data 
sets. We confirmed that the bump feature cannot be de- 
tected today at any significance. We have also checked 
if we can accurately recover w{z) if ACDM was the true 
model for the universe. In this case all the wavelet coef- 
ficients should be zero or the difference from zero should 
not be statistically significant. To test that, we consider 
simulated data for ACDM model and followed the same 
procedure of selecting the first few coefficients with the 
highest signal to noise, and then trying all different com- 



binations of the these selected wavelets. Our MCMC re- 
sults in this case show no evidence of detection of de- 
parture from zero for any of the wavelet coefficients, and 
ACDM is correctly recovered as the best model. 



IV. SEARCHING FOR FEATURES IN 
CURRENT DATA 

In this section, we apply our method to the current 
SNe, CMB and BAO data to search for the deviation of 
■w{z) from —1. We use the recently released "Constitu- 
tion" SNe data sample with 397 data points [25], apply 
the distance priors measured by the WMAP5 team [3], 
and assume a flat universe. 

Table IV shows the main results. We find that using 16 
D4 or Haar wavelet coefficients reduces the and gives a 
hint that w may have deviations from ACDM. However, 
a model with 16 parameters is not preferred by any of the 
information criteria. With D4 wavelets, parameters PI, 
P4, and P8 give the optimal model with a Ax^/A(DoF) 
of —1.6. With Haar wavelets, this number is —1.9 with 
the three wavelet parameters being P2, P4 and P5. We 
also find that none of the wavelet fits are favored by the 
AIC and BIC. 

Plots of w{z) reconstructed using the wavelet method 
with D4 and Haar are shown in Fig. 6. We also plot 
the results obtained using the binning method and the 
CPL parametrization taken from [33]. As one can see, the 
best fit wavelet models show weak hints of dark energy 
dynamics. In particular, w{z) < —1 for z <G [0.25,0.5) 
and w > —1 for z G [0.5, 1], implying that w crosses — 1 
during its evolution. This is consistent with the analy- 
sis of the "Constitution" data set in [33], which reported 
weak evidence for a similar evolution history. Regard- 
less of whether such evolution is real or not, the CPL 
parametrization would not be able to see it. 

A significant deviation from for any of the wavelet 
parameters shows deviation from ACDM. Table V shows 
the eigen-values of principal components for each model 
together with their errors, separately. There is a la and 
2a detection from D4 and Haar respectively, which can 
be interpreted as weak hints that w may have a departure 
from —1. 

We note that although Ref. [33] obtained similar re- 
sults, they used a 5 w bin model. They found this op- 
timal 5-bin model after performing a model-selection on 
10 different binning models, which is computationally ex- 
pensive. With the wavelet approach, we recover the same 
feature in w{z) with only 3 wavelet parameters after run- 
ning MCMC just a few times. This is possible because of 
the multi-resolution nature of the wavelets. We should 
also note that systematic effects like Malmquist bias [34] 
can mimic the same behaviour as dynamical w. There- 
fore, we need to be cautious about making any strong 
conclusions about the equation of state parameter w. 
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FIG. 6: Reconstruction of w{z) from current data using D4 (black solid line), Haar (red dashed line) wavelts, binning (error 
bars) and GPL parametrization (shaded regions). The inner blue error bars (inner shaded dark grey band) and outer green 
error bars (outer shaded light grey band) show the 68 and 95% GL. constraints respectively. 



Model AxVA(DoF) AAIG ABIG 



ACDM 











GDM-D4-16PS 


-0.3 


32.7 


91 


CDM-Haar-16Ps 


-0.4 


31.1 


89 


CDM-D4-1,4,8 


-1.6 


0.9 


13 


GDM-Haar-2,4,5 


-1.9 


0.2 


12.3 



TABLE IV: Differences in the model selection criteria for the 
analysis of the current data. 



Model PGl PG2 PG3 

D4-l,4,8 0.25±1.04 0.36±0.39 0.07±0.13 
Haar-2,4,5 -4.03±1.81 -0.26±0.46 0.10±0.15 

TABLE V: Principal Gomponents for wavelet coefficients of 
optimum w models with current cosmological data 

V. SUMMARY 

We investigated the utility of the wavelet expansion 
method for recovering a local redshift-dependence in the 
dark energy equation of state w{z) with future cosmolog- 
ical data. The advantage of wavelets over other methods 
is that there is no need to assume a scale or a position for 
a feature in w{z). With wavelets, one is able to search for 
small scale features, while looking at the large scale pic- 
ture at the same time. Our results show that a bump-like 



feature in w{z), which is not ruled out by today's data, 
can be recovered using wavelets at a high significance 
with future data. 

We have tried relaxing the flatness assumption and 
considered curvature as an additional parameter. In this 
case, we find that the bump feature we put in becomes 
undetectable, even with the wavelet technique, due to 
the degeneracy between w and the curvature ^Ik- There 
is sufficient freedom for fix and to adjust themselves 
to mimic the effect of a featured w{z) in the luminosity 
distance observations. More specifically, we find that a 
model with w = -1.2 , = 0.137, flbh'^ = 0.0227, 

r^A = 0.6, fiif — 0.1 and h — 0.68 has nearly the same 
goodness-of-fit as the bump model. 

Applying the wavelet method to current cosmological 
data we find a hint of departure from ACDM at a 2cr level. 
The same behavior has been found independently using 
the binning method [33] at the same significance level. 
We have shown that such local features can be found 
more efficiently by using wavelets. The feature we found 
in current data will require further investigation as the 
SNe and other cosmological data improves. 
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